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INTRODUCTION 

This  report  describes  a  program  used  to  calculate  overpressures  in  internal 
explosions  of  aluminized  fuels  in  air.  The  program,  formerly  called  INE,  has  been 
slightly  modified  and  is  called  INAL  (internal  explosions  with  aluminum).  An 
annotated  copy  of  the  INAL  program  appears  in  Appendix  A. 

The  calculations  arc  based  on  the  assumptions  that  the  process  as  a  whole  is 
adiabatic,  that  the  products  arc  uniformly  distributed  through  the  volume,  that  chem¬ 
ical  and  thermal  equilibria  prevail,  and  that  the  ideal  gas  law  may  be  used  throughout. 
The  adiabatic  assumption  in  a  constant  volume  system  is  equivalent  to  the  requirement 
that  there  be  no  net  change  in  internal  energy:  that  is,  that  the  internal  energy  of  the 
products  (34  chemical  species  are  considered)  must  be  equal  to  that  of  the  introduced 
fuel  (inasmuch  as  all  the  other  starting  materials  are  elements  in  thei-  thermochemical 
reference  states). 

A  description  of  the  program  can  rather  logically  be  divided  into  two  segments: 
(1)  the  main  program,  in  which  the  general  procedure  is  to  find,  by  trial  and  error,  a 
product  temperature  for  which  the  adiabatic  condition  holds:  and  (2)  the  subroutine, 
called  Alcal,  which  carries  out  the  equilibrium  calculations  and  then  finds  the  internal 
energy  change  corresponding  to  the  equilibrium  set  of  products. 


MAIN  PROGRAM 


INPUT  SECTION 

Fuel,  enter  formula  and  internal  energy  of  formation.  Computes  formula  mass: 
allows  for  zero  C  or  H  or  no  fuel. 

Cone  (accessible  from  line  1780)  resets  flags,  counters,  and  amounts  of  solids 
to  zero.  Enter  concentration  (kg/m3)  of  fuel  and  of  Ai.  Computes  moles  of  fuel  and 
of  atoms  of  Al,  C,  H,  N,  O.  Computes  total  initial  internal  energy. 
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Temp  (for  entering  trial  temperature  manually)  may  be  accessed  by  use  of 
special  function  key  k4.  After  at  least  two  trials.  Temp  may  be  bypassed  and 
interpolation  used  to  find  the  new  temperature. 


COMPUTATION  SECTION 

Calc  calls  up  the  computational  subroutines  Eq,  Alcal ,  and  Ex. 

Eq  is  a  subroutine  of  the  main  program  to  evaluate  equilibrium  constants  of 
formation  of  each  of  the  34  chemical  species.  Kp  is  first  computed  from  the  stored 
parameters  and  then  converted  to  K„  (expressed  in  mole  numbers).  Kp  is  defined  for 
each  chemical  species  as  the  ratio  of  the  activity  of  the  species  to  the  product  of  the 
activities  of  C,  A1  vapor,  H2,  N2,  arid  02,  each  raised  to  the  power  corresponding  to 
the  stoichiometric  content  of  the  element.  In  converting  the  Kp,  the  activities  of 
gaseous  species,  including  the  elements,  are  changed  to  mole  numbers.  For  the 
condensed  species  the  standard  state  remains  the  pure  phase,  so  that  the  activity  is 
unity,  except  for  C.1 

Alcal  is  the  master  subroutine  which  carries  out  equilibrium  and  energy 
calculations.  Results  are  displayed  as  “dU”  (net)  and  “T  high”  or  “T  low”;  dU  =  0  is 
desired  for  convergence.  A  new  temperature  approximation  is  performed  automatically 
by  interpolation,  based  on  former  T  and  dU  values;  or  else  it  is  entered  manually  with 

k4 . 

Ex  prepares  termination  of  calculation.  Calls  up  Mohum ,  finding  the  total 
moles  of  gases. 


OUTPUT 

Prints  results  of  calculations. 

Cy,  the  molar  heat  capacity  at  constant  volume,  is  computed  by  finding  dU  at 
tv.  nearby  temperatures.  An  isentropic  parameter,  lambda,  is  found  by  using  also 
dn/dT  (change  of  mole  number  of  gas  with  temperature). 

Conclusion:  various  options  allowed,  as  for  new  concentartions  or  a  new  fuel. 


5  Nmt  Weapon*  Center.  Adkbatic  Computation  of  Internal  Blast  from  Atumimm-Cased  Charges  in  Air, 
by  R.  A.  Reinhardt  and  A.  K.  MacDonald,  Nava)  Postgraduate  School  China  Lake,  Calif.,  NWC,  January  1982. 
(NWC  IF  62S7,  publication  UNCLASSIFIED.) 
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ALCAL  SUBROUTINE 

A  leal  is  the  major  computational  subroutine  whose  task  is  to  find  the  numbers 
of  moles  of  the  products  present  at  equilibrium  at  the  selected  temperature.  The 
conditions  to  be  satisfied  (other  than  for  the  trivial  case  of  Ar)  are  the  atom  balance 
conditions  for  C.  H,  N.  O.  and  Al  and  the  establishment  of  chemical  equilibrium 
between  each  compound  and  its  component  elements  in  tiieir  reference  states  at  the 
prevailing  temperature. 

The  master  variables  (all  in  mole  numbers)  are  X  =  y/O;.  V  =  ,/iTJ.  Z  =  v/NT. 
Alv  =  Al  metal  vapor.  Acc  =  activity  of  C  (standard  state  graphite).  Of  these.  Y  is 
always  computed  in  closed  form:  from  one  to  four  of  the  remaining  master  variables 
are  found  as  unknown  parameters,  using  the  Newton-Kaphson  method.  The  actual 
number  of  unknowns  is  equal  to  four,  reduced  by  the  number  of  condensed  phases 
present.  Possible  condensed  phases  are:  AUO;  (solid  or  liquid).  Al  (liquid),  AIN 
(solid),  C  (solid).  AI4C3  (solid).  The  presence  or  absence  of  each  condensed  phase  is 
indicated  by  a  flag,  to  be  set  as  described  later.  The  set  of  condensed  phases  is 
referred  to  as  a  regime:  allowance  is  made  in  the  program  for  about  20  different 
regimes. 

Based  on  the  values  of  the  master  variables  and  the-  Kn  of  formation,  the  mole 
number  of  each  species  is  computed.  Then  the  material  balance  in  the  elements  O.  N. 
C,  and  Al  is  written  in  terms  of  these  mole  numbers.  Thus,  there  results  a  set  of  up 
to  four  simultaneous  non-linear  equations.  It  is  this  set  which  is  used  as  the  basis  of 
the  Ncwton-Raphson  scheme  to  find  the  unknown  parameters. 

At  the  conclusion  of  an  iteration,  the  newly  generated  values  of  the  master 
variables  are  used  to  repeat  the  calculations.  In  favorable  situations  each  iteration 
results  in  improvement  (although  temporary  divergence  sometimes  occurs).  Iteration  is 
repeated  until  errors  in  the  material  balances  are  less  than  one  part  in  ten  thousand. 


INITIAL  APPROXIMATION 

To  begin  the  calculation,  an  initial  approximation  of  the  master  variables  is 
made  using  the  subroutine  Approx.  In  this  approximation  an  arbitrary  hierarchy  of 
oxygen  and  nitrogen  uptake  is  assumed.  Oxygen  is  assumed  to  produce,  in  order.  CO. 
ANO  (g),  ANOj  (c).  HiO.  COi.  and  Oj.  When  there  is  insufficient  O  to  convert  all 
Al  to  ANO;.  AIN  is  assumed  p.esent.  If  there  is  not  sufficient  O  to  convert  al!  C  10 
CO,  AI4C3  is  considered  if  the  temperature  is  low  enough.  This  rather  long  subroutine 
has  given  many  more  satisfactory  initial  conditions  than  she  simpler  scheme  described 
in  Reference  1.  If  is  used  once,  or  at  most  twice,  for  each  new  concentration  that  is 
run. 
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REGIMES  IN  ALCAL 

It  is.  necessary  to  assume  which  condensed  phases  arc  present,  then  perform  the 
calculations  previously  outlined,  and  finally  test  for  the  presence  of  the  condensed 
phases.  The  following  criteria  must  be  met:  (1)  the  quantity  of  the  phase  must  be 
positive  and  (2)  the  formation  constant  must  be  satisfied.  The  tests  are  carried  out 
after  the  convergence  of  the  Newton's  method  calculation.  If  the  tests  fail,  a  different 
regime  (set  of  condensed  phases)  is  tried -the  assumptions  now  based  on  the  currently 
computed  sets  of  mole  numbers. 

An  index  Nw  is  used  to  identify  the  regime  being  considered.  Eleven  differ¬ 
ent  values  of  Nw  arc  allowed,  which  account  for  about  20  different  regimes,  since  for 
most  values  of  Nw  graphite  may  be  present  or  absent.  A  key  to  the  Nw  values  is 
given  at  the  end  of  the  program  in  Appendix  A. 


LABEL  HOM 

The  section  of  A  leal  beginning  at  label  Horn  is  the  setup  for  the  Newton’s 
method  calculation.  (HomO  is  merely  an  early  entrance  into  Ham.)  Depending  on  the 
regime  selected,  li  (the  number  of  variables)  and  Nw  (the  index  identifying  the  regime) 
are  evaluated  by  the  subroutines  liset  and  Nwset.  The  regime  is  established  in  a  variety 
of  ways.  Initially,  Approx  gives  the  first  guess.  Afterwards,  the  subroutine  Almest  is 
used  as  criterion  for  all  solids  except  AI4C3,  which  is  looked  for  at  Uom2.  At  Exit 
and  Exit  I  additional  tests  for  alumina  and  graphite  are  made. 

The  variou.  ibdivisions  of  Horn  are  designated  Nwx.  where  x  is  partly  indic¬ 
ative  of  the  value  of  Nw  for  the  regime  considered.  The  correspondence  is  imperfect 
since  NwJ  includes  the  cases  of  Aw  =1.4,  and  5  and  Nw2  includes  Nw  =  2  and  6. 
In  each  case,  initiii  approximations  arc  set  for  each  variable  to  be  solved,  using  latest 
results  of  the  appropriate  variable. 

Subroutine  Newt  is  called  to  carry  out  the  calculations  (except  Nw  =  Q.  for 
which  solution  in  closed  form  is  possible).  Then  either  a  new  value  of  Nw  is  used  (as 
called  for  by  Alntest  and  liset)  or  the  program  switches  to  Hom2.  As  this  point, 
if  AI4C3  had  not  been  presupposed,  but  was  found  present,  it  is  necessary  to  pass 
once  again  through  Approx.  Otherwise  exit  from  the  subroutine  is  prepared. 

EXIT  FROM  ALCAL 

At  Exit  and  Exit  I  it  is  necessary  to  check  for  negative  amounts  of  alumina 
and  graphite.  If  found,  the  appropriate  flags  must  be  reset  and  a  return  to  Ham  is 
required.  If  those  tests  are  satisfactory,  next,  at  label  Energy ,  the  internal  energy, 
for  each  species  is  calculated  making  use  of  the  stored  parameters  and  the  number  of 
moles  of  the  species.  If  T  =  2315  K  (the  melting  point  of  alumina),  the  relative 
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amounts  of  the  two  phases  of  alumina  are  found  from  the  energy  balance.  Return  to 
the  main  program  then  tikes  place. 


SUBROUTINES  IN  ALCAL 

Alntest  sets  flags  for  AIN,  AI2O3,  Lq  (=  liquid  Al),  and  Gr  (graphite)  based  on 
the  present  value  of  the  formation  expression,  as  related  to  the  formation  constant  for 
the  first  three,  and  on  whether  Acc  has  been  found  less  than  unity  (indicating  no 
graphite). 

liset  and  Nwset  evaluate  //  (the  number  of  variables  to  be  solved  for)  and  Nw 
(the  index  for  the  solids  regime),  based  on  the  flags  set  by  Alntest. 

Newt  solves  the  set  of  simultaneous  non-linear  equations  needed  to  find  the 
values  of  the  master  variables,  using  the  Newton-Raphson  method. 

Fx  gives  the  fitting  functions  for  Newt.  The  subdivisions  are  labelled  Fxy 
where  y  is  a  value  appropriate  to  that  of  Nw.  The  subroutine  generates  a  variable 
designated  Fx  (returned  to  Newt),  which  gives  the  fractional  error  in  the  stoichiometry 
for  whichever  element  is  being  considered  at  the  moment.  Since  Fx  calls  on  Spec, 
mole  numbers  of  all  species  are  found  each  time. 

Diff  first  computes  all  the  master  variables  which  were  not  found  in  Fx  and 
then  calls  Spec  to  find  all  the  mole  numbers  of  the  gaseous  species.  Diff  computes 
mole  numbers  of  the  condensed  species;  and  finally  the  errors  in  the  material  balances 
for  O,  N,  Al.  and  C  are  returned  to  Fx. 

Spec  computes  the  mole  numbers  of  all  gaseous  species,  given  the  current 
values  Of  X,  Z,  Ah-,  and  Acc.  Y  »  y/tfc  is  computed  in  closed  form  in  Spec;  it  is 
needed  for  computations  on  the  hydrogen-containing  species.  FGr  each  species  the  moie 
number  is  computed  by  using  appropriate  values  of  Kn,  the  formation  constant-. 


DIAGNOSTIC  SUBROUTINES 

Sunt  I  is  called  at  the  end  of  each  run  to  give  first  a  check  on  the  material 
balance  in  each  clement  and  tncn:  at  Testk,  a  comparison  of  the  computed  amounts 
with  the  equilibrium  con-rants  of  formation  of  each  condensed  phase  and  of  several 
key  gaseous  species. 

Printt  is  called  whenever  special  function  key  kt)  has  been  depressed  once.  At 
each  emergence  from  Newt  the  mole  numbers  of  all  species  are  given  in  a  condensed 
table,  this  is  followed  by  Suml  and  Testk,  which  give  output  as  previously  described. 
These  various  checks  are  especially  valuable  in  troubleshooting. 
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Depressing  special  function  key  kl  will  cause  execution  of  TRACE 
VARIABLES  Yn(*).  Then  during  each  iteration  in  Newt  the  relative  error  functions, 
used  to  test  convergence,  will  be  displayed.  Each  Yn(J)  must  drop  below  0.0001  in 
absolute  value  before  convergence  is  realized.  This  feature  may  be  turned  off  by 
executing  NORMAL. 
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Appendix  A 
PROGRAM  INAL 

Tills  program  is  written  in  Hewlett-Packard's  HP-BASIC,  and  is  intended  for  use 
on  any  of  the  HP  9845  scries  computers,  in  addition  to  the  main  program,  the  data 
file  *'THD”  and  the  key  file  “INEK"  are  needed  for  operation. 

KEY  FILE  INK 

KEY  • 

Check>-<Cheek«-o> 

-Execute 

*EY  tllACE  VARIABLES  Yn<«> 

-Execute 

KE¥  ?OHI  Ex 
-Execute 

rEY  ^E-STCWE  - 
-Left  arrew 
-Insert  character 

KEY  4 

CONT  Tenp 
-Execute 

KEY  S 

TRACI  PAUSE  Iisct 
-Execute 

KEY  b 

TRACE  PAUSE  S 
-Execute 

KEY  7 

7  .S,3,6 
-Cuntlnue 

KEY  f 

COAT  reel 
-Execute 

KEY  9 

-Clear  line 
LOAD 

KCY10 

■Clear  lire 
SAVE 

KEY11 

■Clear  line 
STORE 

KEY12 

-Clear  line 
EDIT 

KEYlo 

-Clear  line 
EDIT  LINE 

KEV14 

-Clear  line 
LIST 

KEVSS 

-Clear  line 
SCRATCH 

KEYlfc-Cndef  med 
KEY17-Un  defined 
KEY  18-Undcf  med 
KEY19-Und« fined 
KLY2I-Undef ined 
KEY21 -Undefined 
KEY22-Undcf ined 
KEY23 

-117S0 

-Cuntinue 

SFySIWISIS 

KCY£6-Unde fined 
KEY27-Undef ined 

KEY3t-Undefined 
KEY31 -Undefined 
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NWC  TP  6449 
DATA  FILE  THD 

1.  Equilibrium  Constant  of  Formation  Data 


1 

( 

( 

( 

(  1 

-2.M44231721S 

3*95.185437*9 

-61.443183**74 

-7.S97135794(*E-*S 

2.49*78U243»E-*3  -.5 

-4.827587249*5 

24914.8492852 

1.9441*144793 

9.7232S196131E-46 

9.S27274*5371E-»5  -1 

-*.582*6S799Z7 

44171.5374443 

-8.774589**461 

-1 .749738S4821E-M 

4.9892*2484*6E-83  -.5 

2.2S544S587S7 

-E795. 8497374 

27.3534182517 

9.8*7*3939722E-*5 

2.3749393B912E-63  1 

1 

( 

( 

( 

(  ( 

4.59311*31464 

41*8.443*4152 

29.7(8241155 

-1.72129S93489E-I5 

9.17169»97642E-*5  .5 

8. 74863151 113E-42 

2*441.444597 

1.18(739292(1 

-2 . 4347*722*9*E-»S 

6.32488877423E-IS  ( 

3.13219171797 

-12(14.3495425 

19.1339417444 

1.48157371394E-*S 

8.41355776673E-I5  .5 

.749915272297 

•1948.952*9581 

-14.343(133492 

-1.47482381947E-*5 

4.4l2639?*863E-*S  1 

1 

( 

( 

1 

(  « 

-3.»S4*374*471 

133*5.1921*91 

13.8679741254 

-5.19832483354E-I6 

4.7S*84i33874E-*5  -.5 

•2.94319744443 

12919.3248154 

5.24338*84125 

4.78638394752E-I5 

1.8S11593*S84£-*4  -.V 

-4.1476444SS46 

24185.834292 

-7.41178328134 

1 . I4832484614E-I5 

1. 16228* l78*3E-*4  -1 

.71I1I319S414 

-4811.145*3459 

7.783442722*3 

-9 .  M448254772E-46 

3.52l2*847l5*E-(5  1 

1 

1 

( 

1 

i  1 

3.49713844375 

-13438.7244*53 

9.71*861**778 

6.1732276846(E-I6 

fc.933SS5S7*51E-*S  .5 

1 

1 

1 

* 

i  1 

-5.95275584*25 

15458.2581891 

-2*. 84731(4972 

9.l42*6ilt2(6E-aS 

1.71621469496E-»4  -1 

-11 .936351*792 

31(54.2461256 

-37.261473539 

-1.1731*8»8S42E-*4 

2.751863»2573E-*3  -.5 

-41.4169711932 

154392.8*4779 

282.327462644 

1.92799147*57E-»3 

6.138S2»65782E-*3  -3.fi 

-29.4359333294 

1188*7.365935 

-6.1*489862164 

5.*647333«935E-»4 

3.13li(t(8764E-l4  -3.5 

1 

1 

( 

1 

1  1 

.825767716985 

4.76832635592 

-2418.34572*94 

4.65494992481E-85 

1.61847243788E-*2  -.S 

5.19215*83192 

-22218.(4^815 

-11.9581155333 

-1. 4433587 I3*9E-*S 

2.*7761783lilE-l4  .5 

6.52938298135 

-24268.9161784 

-12.471115372? 

-4.297I3319993E-*S 

1 ,86617359728£-*4  .S 

7.16*2273(34 

-29343.3317977 

-.37*425411396 

-6.41291975814E-»5 

9.S477i»77177E-*5  .5 

-24.9*26257877 

71142.893*524 

-49.13*211*137 

-S.19454432333E-*4 

6.44687977746E-*3  -4 

1.61731147515 

-68*6.78*52482 

-12.(46**91672 

-7 . 69 87594781 *E-*6 

S.539*658»187E-*5  1 

-25.74*37539*1 

94491.8463919 

2562.82767313 

2.861852176l4E-(3 

1.15261I13721E-I2  -.S 

2.2214151844 

12*9.(1995973 

219.33532S489 

-S.9968387i43*E-*5 

4.36611S16989E-I5  1 

-2.(3598836359 

6769.63261573 

18.2858441719 

-1.6215553iS26E-*5 

3.6624576232»E-*5  -.5 

2.45411467687 

-11332.242(3*8 

-1S.SS7244H76 

l.3472*7*9929E-*7 

6.45226543872E-K  1 

11.8495*87765 

-41395.(614497 

-16.(531546341 

-1.8S667S2S96*E-*4 

3.897811K914E-I4  1 

26 


NWC  TP  6449 

2.  Internal  Energy  of  Formation  Data 
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3.  Chemical  Species  Symbols 
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A1203ts) 

HP 

CM 

(21 

st 

NCN 

HNCO 
HCO 


C3 
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